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The deconfinement phase transition in Yang-Mills theory with general Lie 
group G 
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We present numerical results for the deconfinement phase transition in Sp{2) and Sp(3) Yang-Mills theories in 
(2 + l)-D and (3 + l)-D. We then make a conjecture on the order of this phase transition in Yang-Mills theories with 
general Lie groups G = SU(N), SO{N), Sp(N) and with exceptional groups G = G(2), F(4), E(6),E{7), E(8). 



1. Introduction and Overview 

Let us consider Yang-Mills theory with gauge 
symmetry group G on a periodic lattice in (d+ 1)- 
D with the Wilson action. This theory is invari- 
ant under the multiplication by an element of the 
center subgroup C(G) of G of all time-like links 
in a given time-slice. This global center symme- 
try is unbroken at low temperatures and it gets 
spontaneously broken at the deconfinement phase 
transition. The Polyakov loop 4> transforms non- 
trivially under this symmetry and thus is an order 
parameter for the deconfinement transition. 

Integrating out the spatial degrees of freedom 
of the Yang-Mills theory, one can write down an 
effective action for $. It describes a scalar model 
with global symmetry C(G) in d-D: the C(G)- 
symmetric confined phase of the gauge theory 
corresponds to the disordered phase of the scalar 
model, while the C(G)-broken deconfined phase 
has its counterpart in the ordered phase. Svetit- 
sky and Yaffe argued that the interactions in 
the effective description are short ranged. Hence, 
if the deconfinement phase transition in the Yang- 
Mills theory is second order, approaching critical- 
ity, the details of the complicated short-ranged ef- 
fective scalar action become progressively unim- 
portant. Only the dimensionality d and the sym- 
metry C(G) of the scalar model are relevant. 



'Speaker at the Conference. 
tOn leave from MIT. 

Work supported by the DOE under the grant DOE-FG03- 
97ER40546 and by the Schweizerischer Nationalfond. 



Thus, one can exploit the universality of the crit- 
ical behavior to use a simple scalar model to ob- 
tain information about the much more compli- 
cated Yang-Mills theory. If the phase transition 
is first order the correlation length does not di- 
verge: there are no universal features and the G- 
symmetric Yang-Mills theory in (d + 1)-D does 
not share the critical behaviour with a C(G)- 
symmetric d-D scalar model. 

Yang-Mills theory on the lattice is naturally 
formulated in terms of group elements while in 
the continuum the fundamental field is the gauge 
potential, living in the algebra. An algebra can 
generate different groups, however it is natural 
to expect that lattice Yang-Mills theories whose 
gauge groups correspond to the same algebra have 
the same continuum limit 2 . Hence, instead of 
SO(N) we consider its covering group Spin(N). 
Keeping this in mind, we look at the center sub- 
groups C(G) of the various simple Lie groups G 



C(SU(N)) = Z(N); 



C(SO{N))- 



C(Sp(N)) 

fz(2); 

>C(Spin(N)) =U{2f- 
U(4); 

C(G(2)) = C(F(4)) = C(E(8)) = {1} 
C(E(6)) = Z(3); C(E(7)) = Z(2) 



= Z(2) (1) 
N odd 
N = 4fc (2) 
A=4fc+2 



(3) 
(4) 



Many numerical simulations in (2 + 1)-D and 
(3 + 1)-D have been performed for SU(N) Yang- 
Mills theory in order to investigate the order of 
the deconfinement transition and - if it is second 
order - to check the validity of the Svetitsky- Yaffe 
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conjecture. A recent study on this subject can be 
found in [3]. The currently known results are: 

• (3 + 1)-D: for TV = 2 - consistent with the 
Svetitsky-Yaffe conjecture - the deconfinement 
transition is second order, in the universality class 
of the 3-D Ising model; for N = 3,4,6,8 the phase 
transition is first order with no universal features. 

• (2 + 1)-D: fluctuations are stronger and the de- 
confinement phase transition turns out to be sec- 
ond order for AT = 2, 3, 4. In agreement with the 
Svetitsky-Yaffe conjecture, the 2-D universality 
classes are, respectively, those of the Ising, 3-state 
Potts, and Ashkin- Teller models. 

The SU(N) branch is not a good choice to 
study the relation between the order of the decon- 
finement transition and the size of the group. In 
fact, when increasing the size N 2 —l, also the cen- 
ter Z(iV) changes. In order to disentangle these 
two features we have considered Sp(N) Yang- 
Mills theory. In this case, the available univer- 
sality class is fixed and the relevance of the size 
of the group on the order of the deconfinement 
transition can be directly addressed. 

2. Sp{N) Yang-Mills Theory 

The matrices U £ SU(2N) with the property 

U* = JUJ + , where J = ia- 2 ®1 (5) 

form the symplectic group Sp(N); o~i is the imag- 
inary Pauli matrix and lis the N x N unit matrix. 
From the definition (|5jl, it follows that the Her- 
mitean generators H of Sp(N) satisfy the con- 
straint H* = -JHJ + = JHJ. Then we can 
write the following general forms for U and H 



U = 



W X 

-X* w* 



H = 



A B 

B* -A* 



(6) 



where A, B,C,D are NxN complex matrices sat- 
isfying: A — A + , B = B T , WW++XX+ =land 
WX T = XW T . Counting the number of degrees 
of freedom, it follows that Sp(N) has N(2N + 1) 
generators and its rank is N. For the center ele- 
ments - since they are multiples of the unit ma- 
trix - it holds that X = and W = W*: hence 
C(Sp(N)) = Z(2). We note that © implies that 
Sp(N) is a pseudo-real group with the special 
cases Sp(l) = SU(2) and Sp{2) = Spin(5). The 



Sp(N) Yang-Mills theory on the lattice can be 
formulated in the usual way in terms of group- 
valued link variables. We have carried out our 
simulations using the Wilson plaquette action. 

3. Numerical Results 

In the next two subsections we report on the re- 
sults of numerical simulations in Sp(2) and Sp(3) 
Yang-Mills theories in (2 + 1)-D and (3 + 1)-D. 

3.1. Sp(2) Yang-Mills Theory 

As a first step we have scanned the expectation 
value of the plaquette from the strong to the weak 
coupling regime. We find no bulk phase transition 
that might interfere with the study of the decon- 
finement transition. In (2 + 1)-D we observe a 
second order deconfinement transition, signalled 
by the broadening of the probability distribution 
of <f> and, hence, by the increase of the Polyakov 
loop susceptibility \ a t criticality. A finite size 
scaling analysis confirms the expectation that the 
universality class is that of the 2-D Ising model. 
Fig. shows the collapse of \ data - collected 
on lattices of different sizes L 2 x 2 and at various 
gauge couplings (3 - on a single curve. We also 
plot rescaled x data for SU{2) Yang-Mills theory 
in (2 + 1)-D: it again deconfines with a second or- 
der transition in the 2-D Ising universality class. 
The two sets agree excellently. In (3 + 1)-D - 
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Figure 1. Finite size scaling function of \j v an d 
7 are 2-D Ising critical exponents; x = (j3/{3 c — 1). 

contrary to what one might have expected - the 
probability distribution of $ in the critical region 
shows the coexistence of the symmetric and of the 
broken phases, indicating that the deconfinement 
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transition is first order. Fig. [21 indeed shows that 
X scales with the spatial volume L 3 at criticality. 
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Figure 2. Scaling of x in Sp{2) Yang-Mills theory 
on L 3 x 2 lattices . We estimate [3 C = 6.4643(3). 

3.2. Sp(3) Yang-Mills Theory 

The results of Sp(2) Yang-Mills theory show 
that in (2 + 1)-D fluctuations are stronger than in 
(3 + 1)-D and the deconfinement transition is sec- 
ond order. Expecting that the larger the group 
the weaker the fluctuations, we have also in- 
vestigated the deconfinement transition in Sp(3) 
Yang-Mills theory. Consistent with this picture, 
we find that (2+l)-D 5p(3) Yang-Mills theory has 
a first order deconfinement transition. The prob- 
ability distribution of $ in the critical region, in- 
deed displays the coexistence of the broken and of 
the symmetric phases. In (3+1 )-D - similar to the 
Sp(2) case - Sp(3) Yang-Mills theory deconfines 
with a first order transition. Fig.[3]shows tunnel- 
ing events between the coexisting symmetric and 
broken phases as a function of Monte Carlo time 
t MC . Sp(3) Yang-Mills theory also has no bulk 
phase transition in (2 + I)-D and (3 + 1)-D. 

4. Conjecture 

Our numerical results show that (2+l)-D Sp(2) 
Yang-Mills theory has a second order deconfine- 
ment transition in the 2-D Ising universality class. 
However, (3 + l)-D Sp(2), (2 + f)-D and (3 + l)-D 
Sp(3) Yang-Mills theories deconfine with a non- 
universal first order transition. Hence, despite the 
fact that a universality class is available, Yang- 
Mills theory can have a non-universal first or- 
der deconfinement transition. A non-trivial cen- 
ter plays no role in determining the order of this 
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Figure 3. (3 + 1)-D Sp(3): tunneling events be- 
tween the symmetric and the broken phases. 

transition. Instead our Sp(N) and the SU(N) |3| 
results indicate that the order of the deconfine- 
ment transition is a dynamical issue related to the 
size of the gauge group. We conjecture that the 
difference in the number of the relevant degrees 
of freedom between the confined phase (color sin- 
glet glueballs) and the deconfined phase (gluon 
plasma) determines the order of the deconfine- 
ment transition. Thus, we expect that in (3+ 1)- 
D only SU (2) Yang-Mills theory has a second or- 
der deconfinement transition; in (2 + 1)-D, due to 
stronger fluctuations, only SU(N), N — 2,3,4, 
and Sp(2) Yang-Mills theories should have sec- 
ond order transitions. According to this picture, 
E(6) and E(8) Yang-Mills theories should also 
have a first order transition due to the large size 
of the groups: 78 and 133 generators, respectively. 
For Yang-Mills theories with trivial center gauge 
groups G(2), F(4), E(8) there is no compelling 
reason for a finite temperature phase transition 
0|. Then, although a first order transition can 
not be ruled out, we expect a crossover. 
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